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SECTION  I 

INTRODUCTION  AND  BACKGROUND 


A.  FLUID  PROPERTIES  AND  THE  EHD  LUBRICATION  PROBLEM 

In  the  elastohydrodynamic  (EHD)  lubrication  problem  the  fluid  undergoes  a 
complicated  stress/deformation  history:  the  stresses  are  transient,  and  have 

both  isotropic  (pressure)  and  shear  components.  (See,  for  example,  reviews  by 
Dowson  and  Higginson  [1]  and  Tallian  [2].)  The  resultant  compressional  and 
shearing  deformations  must  depend  on  the  viscoelastic  properties  of  the  fluid, 
but  in  a rather  complicated,  and  unknown  way.  One  could  pose  the  question: 

Given  a particular  stress  history,  what  are  the  physical  properties  of  the  fluid 
which  must  be  known  in  order  to  determine  the  resulting  deformation  history  (or 
vice  versa)? 

Figure  1 illustrates  the  approach  taken  by  Trachman  [4].  The  pressure 
distribution  (peak  pressure  p0 ) and  geometry  in  the  Hertzian  contact  zone  (of 
width  b)  is  indicated  in  this  figure;  the  fluid  undergoes  a time  varying  pressure 
history  as  it  flows  through  this  contact  region.  By  assuming  a simple  visco- 
elastic model  for  the  liquid,  and  a specific  pressure  dependent  viscosity  func- 
tion, the  "effective"  or  "time  dependent"  viscosity  for  the  liquid  in  the  contact 
zone  is  calculated  and  shown  in  Figure  lb.  By  comparison,  one  sees  that  the 
"effective  viscosity",  hence  shear  traction,  is  much  less  than  by  assuming  an 
equilibrium  fluid  with  no  relaxational  part  to  the  viscosity.  Thus  while  one 
expects  to  enhance  the  viscosity,  the  transient  contribution  is  less  than  the 
steady  state  value,  therefore  partially  compensates.  For  such  calculations  how- 
ever, a more  complete  knowledge  of  the  shear  and  volume  relaxational  behavior 
of  the  liquids  as  well  as  appropriate  time-pressure  equivalence  principles  are 
required.  One  would  expect  that  two  liquids  of  different  chemical  character, 
although  of  the  same  atmospheric  pressure  steady  state  viscosity,  would  exhibit 
a dramatically  different  response  in  the  contact  zone.  Thus,  for  practical  as 
well  as  theoretical  importance,  a fuller  understanding  of  the  viscoelastic 
response  is  required. 


1 


B.  VISCOELASTIC  PROPERTIES  IN  THERMODYNAMIC  EQUILIBRIUM-PHENOMENOLOGICAL 
AND  MOLECULAR  DESCRIPTION 

1.  Behavior  of  the  Linear  Viscoelastic  Parameters 
Tl(uu),  G'(a>),  K'(cu),  k'(u>),  M'  (id)  , M*(tu),  t,,  tv 

In  order  to  phenomenologically  fully  describe  the  mechanical  response  of  a 
liquid  in  thermodynamic  equilibrium  at  a pressure  and  temperature  P,T  we  need  the 
shear  viscosity  function  T)s  (cu) , the  shear  modulus  function  G'liy),  the  volume  vis- 
cosity function  T)v  (uo)  [or  equivalent  K."  (u>)  since  u)T)v  = k"  (cu)  ] , the  bulk  modulus 
function  K'(uj).  It  is  sometimes  convenient  to  approximate  these  functions  by  their 
limiting  values,  and  refer  to  a shear  relaxation  time  ts  and  volume  relaxation  time 
tv  (t  = l/uu)  at  which  a transition  occurs  between  limiting  values. 

To  develop  a physical  feeling  for  the  useful  parameters  in  describing  material 
properties,  it  is  convenient  to  collectively  examine  Figures  3 through  6.  In  Figure  3 
(see  Plazek  [7])  the  specific  volume  (reciprocal  density)  of  the  material  is  shown 
to  vary  with  temperature;  we  shall  allude  to  analogous  variation  with  pressure.  We 
ignore  the  possibility  of  crystallization.  If  we  are  in  the  "liquid"  regime  of  this 
curve,  and  decrease  the  temperature  slightly,  then  one  expects  the  volume  to  decrease 
to  its  new  equilibrium  value  "instantaneously";  certainly  faster  than  our  "time  scale 
of  measurement"  t,  for  which  we  will  assume  r = 1 second.  Thus  one  could  say  that 
the  structural  relaxation  t„,  is  « t . The  same  would  be  true  for  an  increment  in 
pressure.  By  contrast,  in  the  "glassy"  regime  of  this  curve,  in  response  to  a decre- 
ment in  temperature  (or  increase  in  pressure)  the  volume  would  change  very  slowly. 

Thus  the  structural  relaxation  time  tv , the  time  for  the  sample  to  reach  a new  equi- 
librium volume  is  » t ; indeed  tv  could  be  years.  Hence  notice  that  the  glass  tem- 

perature is  not  unique  but  depends  on  whether  tv  » t or  tv  « t,  namely  depends  on 
our  choice  for  the  time  scale  of  the  measurement. 

Consider  as  an  illustrative  example  of  pressure  dependence  of  the  bulk  dynamic 
compliance,  B*(B*  = 1/K*,  where  K*  is  the  bulk  modulus),  the  reduced  variable  curve 
of  Marvin  [12],  Figure  5b.  The  measurement  time  scale  is  taken  as  t ~ 1/1000  sec; 

at  a reduced  pressure  P,  =-2000  bars  (or  at  high  temperatures),  tv  « t and  we  obtain 

the  liquid-like  compliance  B0 , whereas  at  a reduced  pressure  Pr  = +2000  bars,  tv  » t 
and  then  we  obtain  the  glass-like  compliance  B^.  In  these  examples,  the  t of  the 
m isurement  is  kept  fixed  and  t„  is  changed  by  pressure  or  temperature'.  Alternatively, 
if  the  P and  T are  held  constant  as  in  Figure  5a  and  only  i is  changed,  then  one 
measures  a liquid-like  compliance  b,  at  t » tv , and  glass-like  Be  for  i « tv  . This 
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figure  illustrates  that  tv  can  be  determined  at  each  fixed  T and  P if  t is  variable 
over  a wide  time  (frequency)  range.  We  now  conclude  that  even  the  "liquid-like" 
regime  of  Figure  3 can  behave  "glass-like"  if  t is  short  enough. 

Hence  the  viscoelastic  properties  of  the  lubricant  Be(P),  Bg (P)  and  tv (P) 
must  be  determined  as  a function  of  pressure.  And  more1. 

Turning  next  to  the  viscoelastic  response  of  a lubricant  to  (dynamic)  shear 
deformation,  Figures  4a  and  b indicate  additional  parameters  required  to  describe 
the  fluid.  There  are  two  characteristic  relaxation  times+  tc  and  ts  if  the 
lubricant  contains  polymers,  through  molecular  models  these  are  associated  with 
conformational  change  of  the  polymers  tc,  and  orientational  changes  of  the  solvent 

molecules  ts  . Thus,  for  t » tc  > ts  (or  « u>0  < u)s  ) , the  lubricant  behaves 

liquid-like  with  a shear  viscosity  and  "zero"  shear  modulus;  for  t0  > t > t6 , 

the  viscosity  drops  to  a new  level  of  T],*,  and  the  shear  modulus  G'  increases  to 

rubber-like  values;  while  for  t < ts  the  viscosity  drops  further  and  G'  reaches 
a glass-like  value. 

Hence  the  shear  viscoelastic  properties  of  the  lubricant  Tb(P),  Tj^P),  Gro(P), 
tc  and  ts  must  also  be  determined  as  a function  of  pressure. 

As  an  aside,  it  should  be  mentioned  that  a polymeric  fluid  being  sheared 
will  also  generate  normal  stresses  which  tend  to  push  the  contact  surfaces  apart. 
The  magnitude  of  these  normal  stresses  are  related  to  and  rc . 

2.  The  Free-Volume  Picture:  Concepts  of  Time-Temperature  Superposition 

One  unifying  model  of  liquids  which  has  had  a large  degree  of  success  in 
tying  together  temperature  and  pressure  dependent  properties  has  been  the  free 
volume  model  of  Doolittle  [8].  The  idea  is  that  the  actual  volume  occupied  by 
the  molecules  is  less  than  the  actual  sample  volume  v,  so  that  the  difference 
is  the  free  volume  vf . The  larger  the  fractional  free  volume  f = vf /v,  the 
easier  it  is  for  molecules  to  slide  past  each  other  or  to  rearrange  conformations. 
Thus  Doolittle  postulated  that  the  temperature  and  pressure  dependence  of  f 
would  be  linear: 

f = f0  + ar  (T  -Tq)  + gr  (P  -P0)  (1) 


In  every  case,  by  "characteristic  relaxation  time"  we  mean  the  mean  value  of 
the  relaxation  time  distribution  function. 
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and,  further:  that  at  two  different  temperatures  T and  T0 , each  characteristic 
relaxation  time  would  change  in  the  same  ratio  as  the  zero  shear  viscosity, 
then  a temperature  shift  factor  is  defined: 

_ Vt 


= t8/t 

t=/t„ 


Vt0 


from  which  he  obtains 


a'  = K?  ' -t) 


This  equation  was  extended  by  Ferry  [15]  to  the  WLF  equation 

9m  .T  - 


C2  +T  -T0 


(2) 


(3) 


(4) 


which,  in  combination  with  Eq.(2)  has  formed  a time-temperature  superposition 
principle  that  has  been  successful  when  applied  to  conformational  relaxation 
times.  The  inference  of  such  an  equation  is  that  the  character  of  the  relaxation 
mechanism  does  not  change  with  temperature  (pressure),  but  only  the  magnitude  and 
characteristic  time  scale  will  change.  This  is  not  true  for  all  relaxation  mech- 
anisms. The  utility  of  a time-temperature  superposition  principle  is  that  pre- 
dictions can  be  extrapolated  to  time  (or  temperature)  ranges  other  than  those 
measured . 

3.  The  Free-Volume  Picture:  Concepts  of  Time-Pressure  Superposition 

Time-pressure  superposition  is  far  less  developed  for  large  pressure  ranges. 
Ferry  [15]  has  extended  the  concepts  above  to  obtain: 

_ (B/2.303  f0)(P  -P0) 


log  ap  = 


f0/3f  - 0*  *po 


(5) 


The  free-volume  model  would  say  that  only  the  final  volume  is  important  regard- 
less of  the  temperature  or  pressure  history  route.  This  idea  needs  further 
testing.  Other  approaches  have  been  developed  by  Cohen  [9]  and  by  Gibbs  [10]. 

C.  VISCOELASTIC  RESPONSE  TO  THE  TRANSIENT  (NONEQUILIBRIUM)  EHD  SITUATION 

This  subject  of  the  relation  of  equilibrium  properties  to  transient  response 
has  not  thus  far  been  studied  extensively,  and  therefore  is  part  of  our  experi- 
mental goals.  One  useful  view  taken  by  Harrison  [11a]  and  by  Trachman  [4]  uti- 
lizes free  volume  concepts  discussed  above.  The  central  argument  can  be  viewed 
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as  if  volume  is  the  unifying  concept,  then  pressure  (and  temperature)  history 
may  be  put  aside,  and  the  response  may  be  determined  by  integration  over  volume 
history,  with  the  assumption  of  equilibrium  properties  at  each  segment  of  the 
history. 

D.  EXPERIMENTAL  TECHNIQUES 

1.  Mechanical  Spectroscopy:  Steady  State,  Low  Frequency,  Ultrasonic 

The  mechanical  relaxation  spectra  of  interest  for  lubrication  liquids  can 
reach  from  steady  state  (D.C.)  to  lO10  Hz.  The  part  of  the  spectrum  of  interest 
will  of  course  change  for  each  pressure  (and  temperature).  As  indicated  earlier 
the  functions  which  must  be  known  are  the  shear  viscosity  and  modulus,  and  the 
bulk  viscosity  and  modulus.  While  a variety  of  mechanical  techniques  extending 
from  1 Hz  to  10^  Hz  have  been  used  in  this  laboratory  (Figure  2)  for  liquids  at 
atmospheric  pressure,  not  all  of  them  are  directly  applicable  to  high  pressure 
techniques.  Some  exceptions  are  applications  of  similar  methods  to  high  pressure 
problems  by  Barlow  [11],  Marvin  [12]  and  by  Lamb  [13]. 

The  primary  approach  taken  in  our  studies,  is  to  construct  a practical  high 
pressure  cell  optimized  so  that  the  widest  frequency  range  of  measurements  can 
be  made  in  the  same  cell  - both  under  equilibrium  and  transient  pressures.  Be- 
cause of  the  limit  of  small  sample  volume,  we  emphasize  the  techniques  which 
utilize  light  scattering,  and  those  which  couple  mechanical  viscoelastic  techniques 
with  light  scattering  - rheo-optic  techniques. 

2.  Light  Scattering  - Twelve  Years  of  Progress 

With  the  advent  of  lasers,  fast  photomultipliers,  highest  quality  Fabry- 
Perot  interferometers,  high  speed  (computerized)  autocorrelators  - material 
property  information  - extracted  by  light  scattering  has  developed  very  far 
indeed  these  past  twelve  years.  Here  we  will  only  list  a few  techniques  relevant 
to  lubrication,  a few  developed  in  our  laboratory  for  the  past  eight  years,  many 
developed  elsewhere. 

The  general  features  of  the  light  scattering  experiment  is  found  in 
Figure  8.  Light  from  a laser  is  focussed  into  a narrow  beam  (circa  0.1  mm)  at 
the  liquid  sample,  and  the  light  scattered  at  any  set  scattering  angle  9 is 
collected  by  a photomultiplier.  The  "liquid"  scattering  cell  may  be  one  of  the 
new  high  pressure  cells  of  Figures  9 or  11,  or  our  older  high  pressure  cell  (22) 


of  Appendix  1.  The  scanning  photometer  (26)  paper  provides  a more  detailed 
description.  Optical  detection  can  be  direct  as  shown  in  Figure  8 for  autocor- 
relation and  cross-correlation  function  analysis  of  scattered  intensity;  or  a 
Fabry-Perot  Interferometer  (FP)  may  be  inserted  between  liquid  and  photomulti- 
plier (PM)  for  Brillouin  peak  and  Rayleigh  Wing  analysis;  filters  for  fluores- 
cence studies  and  analyzers  for  depolarized  scattering  are  also  introduced  be- 
tween liquid  and  PM. 

(i)  Rayleigh  line.  Figure  7 shows  a light  scattering  spectrum  from  Benzene 
at  9 = 90°  taken  by  Cummins  [16]  using  a FP.  The  central  line  is  called  the 
Rayleigh  line  because  the  center  frequency  is  the  same  as  that  of  the  laser. 
However,  whereas  the  laser  line  width  is  a few  Hz,  this  Rayleigh  line  is  ~ 10s  Hz 
and  is  due  to  nonpropagating  entropy  fluctuations,  i.e.,  changes  in  local  order- 
ing of  the  molecules  or  rotational  ordering  of  the  molecules  (when  measured  in 
HV) . Notice  VV  refers  to  the  polarized  component  of  the  light,  HV  to  the  de- 
polarized component.  Speed  of  molecular  motion  is  related  to  line  width. 

(ii)  Brillouin  lines  W polarization.  The  lines  on  either  side  of  the 
central  line  in  Figure  7 are  due  to  pressure  fluctuations  inherent  in  all  ma- 
terials; this  becomes  equivalent  to  an  ultrasonically  generated  pressure  dis- 
turbance which  of  course  propagates  as  a wave.  Such  waves  are  spontaneously 
present  at  all  frequencies  (and  wavelengths),  so  that  one  selects  the  acoustic 
wavelength  by  setting  9,  the  scattering  angle,  and  measures  the  frequency  shift 
with  respect  to  the  central  line  to  obtain  the  velocity.  Together  with  the 
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density  of  the  fluid,  this  yields  the  longitudinal  modulus  M,  where  M=K  + j G, 
see  for  example  Dill  [23]  or  McSkimin  [24].  A more  detailed  discussion  of  the 
nuances  in  extracting  K and  G,  as  required  for  lubricant  properties  is  omitted 
here.  The  width  of  these  Brillouin  peaks  yields  the  wave  attenuation,  or  dynamic 
viscosity  T]H  . 

(iii)  Brillouin  lines  HV  polarization.  If  the  frequency  and  pressure  are 
such  that  the  material  behaves  glass-like,  and  can  sustain  a shear  stress,  then 
a transverse  velocity  and  attenuation  corresponding  to  Go,  and  Tlo,  can  also  be 
determined . 


(iv)  Quasi-elastic  scattering.  When  polystyrene  spheres  of  size  < 
are  added  in  minute  concentration  to  the  liquid,  they  become  the  primary  scatter- 
ing source,  and  their  relative  motions  give  rise  to  time  varying  interference  of 
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the  light  reaching  the  PM  (no  FP) . The  resulting  photocount  fluctuations  are 
analyzed  by  the  computer/autocorrelator  as  seen  in  Ref. 26  and  Appendix  1.  The 
autocorrelation  function  yields  diffusion  coefficient  of  the  spheres,  hence  one 
obtains  the  viscosity  (low  frequency  limit)  T]0  of  the  liquid. 

(v)  Quasi-elastic  scattering.  A new  technique  under  development  is  to  re- 
place the  polystyrene  spheres  used  above  with  colloidal  magnets  and  force  me- 
chanical oscillation  at  medium  frequencies  (1  Hz  to  10  KHz)  by  an  external  mag- 
netic field.  The  motion  is  analyzed  by  a cross -correlation  technique  related  to 
the  one  used  above. 

E.  OBJECTIVES  OF  THIS  RESEARCH  PROGRAM 

In  elastohydrodynamic  (EHD)  lubrication  the  liquid  lubricant  is  subjected 
to  a complex  transient  stress  history  in  which  the  isotropic  pressure  components 
may  reach  beyond  200,000  psi  in  durations  of  the  order  of  1 millisecond.  One 
may  ask:  How  do  the  viscoelastic  characteristics  of  a fluid  influence  the  stress 

built  up  in  this  lubrication  situation:  Therefore  our  long  term  objectives  are: 

(i)  To  measure  the  rheological  properties  of  a variety  of  liquid  lubricants 
over  the  full  frequency  range  of  viscoelastic  response  from  1 Hz  to  1010  Hz; 
over  a pressure  range  to  300,000  psi;  and  over  a temperature  range  from 
-100° C to  +300°C.  Knowledge  of  the  viscoelastic  response  is  important  both 
under  equilibrium  and  transient  pressure  and  temperature  histories,  as  well 
as  in  uniform  and  nonuniform  pressure/temperature  fields. 

(ii)  Phenomenological  theories  are  to  be  explored  which  can  relate  the  above 
equilibrium  state  viscoelastic  measurements  to  the  stress/strain  response  of 
a liquid  under  transient  elastohydrodynamic  lubrication  conditions.  In 
other  words,  the  studies  should  lead  to  time-temperature  and  time-pressure 
superposition  principles  applicable  under  transient  and  nonuniform  conditions. 

(iii)  Molecular  theories  are  to  be  examined  which  relate  molecular  struc- 
ture to  the  measured  macroscopic  response;  thereby  establishing  possible 
relations  between  lubricant  properties  and  molecular  structure. 

(iv)  Experimental  procedures  in  this  laboratory  are  based  primarily  on 
rheo-optic  techniques;  that  is,  application  of  stresses  by  external  mechani- 
cal or  by  spontaneous  means,  and  measurement  of  response  by  optical  (laser 
light  scattering)  means. 
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(v)  Selection  of  materials:  After  consultation  with  Mr.  F.C.  Brooks  and 

Dr.  J.F.  Dill  it  was  established  that  the  fluid  PR-143  AC  would  be  most 
suitable  for  initial  studies,  and  Emery  2958  and  Stauffer  704  would  be 
appropriate  follow-on  candidate  fluids.  In  this  discussion  it  was  estab- 
lished that  the  most  meaningful  progress  could  be  made  if  coordination 
between  laboratories  was  established  on  the  basis  that  each  measured  some 
identical  fluids . Of  course,  an  understanding  of  the  relationship  between 
chemical  structure  and  mechanical  properties  of  lubricants  remains  a 
continuing  objective. 
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SECTION  II 

EXPERIMENTAL  PROGRESS  AND  RESULTS 

A.  THE  LIGHT  SCATTERING  SPECTROMETER 

Our  light  scattering  spectrometer  is  in  the  process  of  being  upgraded  to 
add  the  capability  of: 

(i)  Brillouin  line  width  measurements  for  polarized  and  depolarized 

components . 

(ii)  Scattering  from  magnetic  spheres  in  linear  motion  in  the  liquid. 

(iii)  Scattering  from  magnetic  spheres  in  oscillatory  motion  in  the 

liquid.  These  are  additional  to  the  capability  hitherto  available  of: 

(iv)  Rayleigh  (quasi-elastic)  scattering. 

(v)  Scattering  from  suspended  polystyrene  spheres  in  Brownian  motion. 

In  addition,  for  measurement  of  pressure,  we  must  add  the  capability 

to  measure: 

(vi)  Fluorescence  spectral  shift  from  ruby  crystals. 

A multiport  microscope  system  has  been  designed  to  accommodate  this  variety 
of  optical  methods  from  the  small  volume  of  the  high  pressure  cells.  The  major 
components  of  this  new  system  have  now  been  procured  or  were  available  from  the 
previous  system.  It  is  useful  to  summarize  the  design  considerations  and 
present  ststus  of  each  of  the  major  components  of  the  light  scattering  spectro- 
meter system. 


1.  The  Light  Scattering  Photometer 

The  restrictions  imposed  by  a small  volume  high  pressure  cell  (see  dis- 
cussion of  cell),  together  with  the  requirement  for  several  experiments  have 
led  us  to  the  design  of  a multiport  microscope  light  scattering  photometer  to 
replace  our  current  apparatus.  The  microscope  will  have  three  ports  as  seen  in 
Figure  12:  one  for  direct  visual  or  video  scan  observation  of  liquid  to  glass/ 
crystalline  transformations  (and  kinetics  thereof)  under  high  pressure;  one  port 
for  quasi-elastic  light  scattering  (Rayleigh  line)  data  analysis  including 
oscillating  magnetic  spheres  and  for  Brillouin  and  Rayleigh  wing  spectral  analy- 
sis; and  one  for  determining  local  pressure  by  the  fluorescent  spectral  shift 
of  ruby  as  measured  with  an  optical  scan  spectrometer. 
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For  the  diamond  anvil  high  pressure  cell,  surface  scattering  can  be  a 
serious  problem.  In  our  previous  experience  with  scattering  from  capillary 
tube  cells,  we  were  able  to  first  keep  the  surface  scattering  localized  by  ad- 
justment of  the  angle  between  the  cell  and  incident  beam,  and  then  adjustment 
of  the  angle  and  volume  of  observation  so  that  the  surface  regions  were  excluded 
from  view.  Since  a (estimated)  + 20°  rotation  of  optical  access  to  the  cell  is 
possible,  we  must  provide  for  these  angular  adjustments  to  minimize  the  amount 
of  surface  (liquid-diamond  interface)  scattered  light  that  enters  the  optical 
system. 

To  accommodate  the  remaining  surface  effects:  (a)  For  autocorrelation  the 

surface  will  act  as  a reference  signal  (i.e.,  no  frequency  shift),  therefore 
the  method  will  become  hetrodyne  instead  of  homodyne;  (b)  for  Brillouin  scatter- 
ing the  unshifted  line  can  be  tuned  (with  our  new  laser  - see  below)  to  an 
Iodine  (I2)  absorption  band.  This  means  the  I2  acts  as  a sharp  band  reject 
filter  of  width  900  MHz;  (c)  for  fluorescence  the  additional  scattering  is  less 
important.  The  air  suspended  granite  optical  table  remains  excellent  for  the 
new  system.  New  mounting  fixtures  will  be  required  for  the  new  laser,  Fabry- 
Perot,  multiport  microscope  and  goniometer. 

2.  The  Argon  Laser 

One  of  the  major  new  components  is  the  Argon  Laser  to  replace  our  present 
HeNe  laser.  For  Brillouin  line  width  measurements,  particularly  at  small  scat- 
tering angles  and  where  correlation  techniques  and  the  Fabry-Perot  interferometer 
techniques  can  overlap,  it  is  essential  that  the  laser  source  be  as  stable  (in 
frequency)  as  possible.  Amplitude  stability  is  required  for  the  quasi-elastic 
(intensity  fluctuation)  experiments.  Towards  this  end,  the  following  test  for 
amplitude  and  frequency  stability  was  devised. 

The  laser  etalon  was  first  adjusted  so  that  the  single  frequency  of  the 
laser  corresponds  to  the  absorption  peak  in  an  Iodine  absorption  cell  and  the 
intensity  Ip  is  monitored  on  an  oscilloscope.  Then  the  etalon  is  then  adjusted 
to  the  shoulder  of  the  Iodine  absorption  and  the  intensity  Ie  is  monitored. 

The  percentage  fluctuations  in  Ip  are  due  to  amplitude  noise,  while  fluctuations 
in  I(  contain  both  amplitude  and  frequency  noise.  The  spacing  between  longi- 
tudinal modes  of  the  laser  were  used  for  frequency  calibration. 
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The  two  commercially  available  argon  lasers  considered  most  stable,  Spectra 
Physics  Model  165  and  Lexel  Model  85,  were  tested  in  this  way.  On  the  basis  of 
performance,  the  Lexel  laser  was  found  preferable. 

A more  powerful  version  of  the  Lexel  85,  the  Lexel  95,  was  finally  selected 
and  has  been  delivered.  We  are  attempting  to  pre-chill  the  laser  cooling  water 
to  achieve  laminar  flow  of  the  coolant,  since  acoustic  noise  generated  by  turbu- 
lent flow  is  partially  responsible  for  the  noise  in  the  laser  frequency  spectrum. 

3.  The  Fabry-Perot  Interferometer 

Another  major  component  is  the  Fabry-Perot  interferometer  for  Brillouin 
shift  and  line  width  measurements.  After  comparison  of  specifications  of  Tropel 
and  Burleigh  units,  the  Burleigh  Fabry-Perot  Interferometer  Model  RC10  with 
multipass  option  was  selected.  Aside  from  stability,  we  concluded  that  the 
multipass  feature  of  this  unit  will  be  very  helpful  in  the  reduction  of  the 
contribution  of  scattering  from  the  surfaces  of  the  diamond  cell.  Agreement 
was  reached  with  RPI  for  the  purchase  of  this  unit,  and  the  unit  has  been 
delivered . 

The  unit  is  being  set  up  and  interfaced  to  the  computer.  The  electronics 
for  individual  adjustment  of  the  three  piezoelectric  drivers  and  the  high  voltage 
sweep  circuit  must  still  be  constructed  and  interfaced  to  the  digital  to  analog 
drives  from  the  computer  (Figure  10).  Also  the  computer  program  for  tuning  the 
interferometer  and  for  automatic  frequency  drift  correction  must  still  be 
written. 

4.  Photodetection 

The  photodetection  system  based  on  our  present  thermoelectrically  cooled 
Amperex  56  TVP  photomultiplier  will  be  adequate  for  the  present;  a lower  back- 
ground ITT  FW130  phototube  may  be  desirable  at  a later  stage.  An  adapter  to 
mechanically  couple  this  phototube  housing  to  the  collimator  collection  optics 
of  the  Fabry-Perot  must  be  constructed. 

5.  Data  Acquisition  System 

The  present  computerized  data  acquisition  and  analysis  (autocorrelation) 
system  as  seen  in  Figure  12,  continues  to  work  well.  Only  more  software  will 
be  required  for  the  new  functions. 
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No  additional  hardware  is  required.  The  same  software  auto-  and  cross  - 
correlation  programs  will  be  used,  and  only  an  automatic  sweep  normalization 
program  has  to  be  written  for  the  new  Fabry-Perot. 

6.  Other  Considerations 

Several  of  the  components  of  our  earlier  light  scattering  system  are  being 
incorporated  in  the  new  Diamond  Cell  Light  Scattering  geometry.  We  are  still 
in  the  process  of  assembly,  so  that  some  minor  components  must  still  be  selected, 
and  several  mounting  fixtures  must  be  constructed  in  the  workshop. 

For  the  high  pressure  dynamic  viscosimeter  method,  magnetic  spheres  are 
suspended  in  the  lubricant  in  the  diamond  cell,  and  these  are  oscillated  by  an 
external  magnetic  field. 

We  have  explored  a variety  of  colloidal  magnetic  particles  and  tiny  iron 
spheres  for  use  in  the  high  pressure  cell  dynamic  viscosity  technique.  The 
colloidal  suspensions  compatible  with  the  lubricants  to  be  studied  have  shown 
a tendency  to  aggregate.  The  tiny  iron  spheres  seem  to  work  well  in  the  lubri- 
cants. These  were  available  from  Whittaker  Corporation  in  250  p.  and  50  p 
diameters,  suitable  for  our  experiments.  They  can  be  individually  handled.  We 
are  testing  a variety  of  magnetic  field  drive  geometries. 

B.  HIGH  PRESSURE  CELLS 

We  have  decided  to  pursue  two  high  pressure  cell  designs  in  this  laboratory, 
the  cross  cell  and  the  diamond  anvil  cell. 


1.  The  Cross  Cell  Design 

Earlier  measurements  [22]  on  light  scattering  from  suspended  spheres  were 
made  in  this  laboratory  (Appendix  A)  over  the  pressure  range  up  to  35,000  psi. 
The  cell  worked  well,  but  the  range  is  less  than  pressures  generated  in  the  EHD 
regime,  so  a new  cross  cell  of  similar  design  was  proposed  to  extend  the  range 
beyond  150,000  psi.  In  this  modified  design  (Figure  9),  the  sample  volume  was 
substantially  reduced  so  that  pressure  transmission  would  not  be  impeded  by  the 
transition  of  the  sample  from  liquid-like  to  glass-like  behavior,  and  so  that 
the  same  valve  stem  would  serve  as  the  compression  piston.  The  latter  would 
then  permit  transient  pressure  pulses  to  be  applied. 
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The  major  advantage  of  this  cross  cell  design  is  that  scattering  can  be 
observed  at  90°  thereby  eliminating  the  surface  scattering  contribution.  The 
disadvantage  is  that  material  failure  limits  the  highest  achievable  pressures. 

2.  Diamond  Anvil  Cell 

The  basic  design  of  the  reho-opcic  light  scattering  high  pressure  cell 
instrument  is  centered  on  microscopic  observation  of  the  fluid  under  high  pres- 
sure. The  use  of  a tiny  volume  as  in  the  diamond  cell  (see  Figure  11)  seems 
essential  if  one  is  to  achieve  pressures  in  the  300,000  psi  range,  since  the 
transition  from  liquid-like  to  glass-like  behavior  of  the  sample  would  resist 
the  pressure  build-up.  Furthermore,  the  cell  volume  becomes  orders  of  magnitude 
closer  to  the  volume  in  EHD  contacts,  so  that  one  brings  the  laboratory  experi- 
ment much  closer  to  the  real  contact  situation.  The  process  of  generating  high 
pressures  (and  impulses)  becomes  simplified,  a single  hand  screw  is  sufficient 
to  generate  the  desired  pressures.  The  additional  complexities  that  arise  from 
the  microscopic  volume  are  connected  with  sample  loading,  measurement  of  pres- 
sure, and  surface  scattering.  Furthermore,  microscope  optics  are  required  for 
the  small  sample  volume,  and  were  described  above. 

A diamond  anvil  cell  designed  by  Dr.  A.  Van  Valkenberg,  while  at  the 
National  Bureau  of  Standards,  was  ordered  and  delivered  from  High  Pressure 
Diamond  Optics,  Incorporated.  Both  optical  and  sample  handling  procedures  are 
currently  being  developed  for  these  samples  of  volume  circa  1CT4  cc. 

Pressure  and  pressure  distribution  measurement  can  be  determined  by  the 
shift  in  the  fluorescence  spectra  of  ruby  crystals  suspended  in  the  fluid  [18. 

19]  and  locally  monitored  by  the  microscope  optics.  Design  is  underway  for  the 
fast  pressure  pulse  triggering  device. 

C.  ULTRASONIC  SHEAR  WAVE  APPARATUS 

1.  Introduction  and  Preliminary  Results 

There  are  several  atmospheric  pressure  reference  state  measurements  that 
are  required  for  proper  interpretation  of  higher  pressure  data.  One  of  the  more 
important  is  the  high  frequency  limiting  shear  viscosity  and  modulus  (transverse 
wave)  and  relaxation  time(s).  Steady-state  shear  viscosity  (low  frequency  limit) 
data  is  available  for  the  test  fluids,  however,  the  directly  measured  high  fre- 
quency limiting  shear  modulus  appears  to  be  lacking.  Since  Brillouin  scattering 
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from  shear  waves  is  difficult  to  observe,  the  high  frequency  limit  and  T]^ 
must  be  calculated  from  the  longitudinal  modulus  with  some  assumptions.  There- 
fore, we  assembled  auxiliary  ultrasonic  delay  line  apparatus  for  direct  measure- 
ment of  G'  and  Tl'  even  though  it  is  only  useful  at  1 atmospheric  pressure. 

This  apparatus  is  currently  being  calibrated  in  the  frequency  range  2 to  50  MHZ 
and  over  a temperature  range  from  -10°C  to  +70°C.  Thus  this  apparatus  serves 
two  purposes:  (i)  to  provide  the  high  frequency  shear  moduli  and  relaxation  times 
mentioned  above  for  bulk  lubricants  and  (ii)  to  provide  the  limiting  shear  re- 
laxational  behavior  for  lubricants  with  polymer  additives  (i.e.,  a measure  of 
the  configurational  relaxation  contribution  of  the  polymer  additive). 

This  apparatus  has  been  constructed,  tested,  calibrated,  and  we  are  now  in 
the  process  of  measuring  Krytox  143. 

Calibration  curves  have  been  obtained  for  several  fluids  which  we  believe 
to  be  Newtonian  up  through  the  frequencies  of  interest. 

Preliminary  results  have  also  been  obtained  on  solutions  of  polystyrene  in 
a di-2-ethyl-hexyl  phthalate  fluid.  These  results  are  relevant  to  the  role  of 
polymer  additives;  and  indicate  that  at  high  frequencies  the  viscosity  of  the 
solution  does  not  become  that  of  the  solvent.  Namely,  theory  has  assumed  that 
when  conformational  relaxation  of  the  polymer  has  been  "frozen  out",  then  only 
the  viscosity  of  the  solvent  would  remain.  This  picture  is  incomplete.  More 
extensive  data  and  analysis  of  this  phenomenon,  called  the  T]oo  problem,  is  in 
progress . 

2.  Description  of  Apparatus 

The  description  of  this  apparatus  may  be  found  in  Appendix  B,  coupled 
with  Figures  10  and  13.  The  derivation  of  basic  device  equations  follow. 


lit. 
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3.  Derivation  of  Basic  Surface  Wave  Loading  Equations 

for  Delay  Line 

In  this  experiment  the  main  interest  is  measuring  the  complex  shear  modulus 
p,*  of  the  liquid  sample  which  is  in  contact  with  major  surfaces  of  the  delay 
line.  The  equations  for  propagation  of  ultrasonic  waves  in  a plate  are  developed 
in  this  section  and  then  restricted  to  the  case  of  interest,  namely,  the  lowest 
order  transverse  mode.  This  mode  is  often  called  the  Meitzler  mode  [1],  We 
obtain  solutions  for  the  velocity  and  attenuation  of  this  wave  in  the  presence 
and  absence  of  liquid  at  the  interface.  Equations  for  the  relationship  between 
velocity  and  attenuation  change  and  the  mechanical  impedance  and  the  shear  modulus 
of  the  liquid  are  also  obtained. 

Consider  an  infinite  plate  (Figure  13)  of  thickness  2b  in  a coordinate  system, 
where  it  is  bounded  by  the  major  planes  at  X = b and  X = -b  with  a liquid  sample. 
When  an  ultrasonic  transverse  wave  propagates  down  the  plate  in  the  z-direction 
a wave  is  also  generated  in  the  liquid  because  of  the  no-slip  condition  between 
solid  plate  and  liquid  at  the  interface. 

The  steady-state  sinusoidal  excitation  with  the  appropriate  boundary  con- 
ditions results  in  a set  of  homogeneous  equations  whose  solutions  have  discrete 
modes  of  wave  propagation. 

For  the  delay  line  used  here,  the  dimensions  and  polarization  of  the  piezo- 
electric (P.Z.T.  ceramic  transducer  material)  are  chosen  so  that  only  the 
principal  shear  mode  (Meitzler  mode)  of  propagation  is  excited  while  higher 
shear  modes  are  suppressed  (when  the  plate  is  in  vacuum).  It  is  very  useful  to 
have  only  the  nondispersive  Meitzler  mode  SS(O)  because  its  velocity  is  inde- 

' pendent  of  frequency  and  the  complication  of  having  several  modes  of  different 

j 

I . velocities  can  be  avoided. 

* 

* 

Even  in  the  presence  of  liquid  on  the  major  surfaces  it  will  be  assumed 
that  only  the  same  SS(O)  mode  is  present  since  the  mechanical  impedance  of  the 
liquid  is  much  less  than  the  solid. 

The  equation  of  wave  propagation  for  an  isotropic  Lomogeneous  linear  visco- 
elastic material  is  generally  written  as  follows: 

liVu  + (X  +M,*)  v(v  • U)  = P — r , (l) 


where 


p*  - complex  shear  modulus 
\ - other  Lame  coefficient 
u - particle  displacement  vector 
p - density. 

If  one  uses  Hamilton's  decomposition  method  by  defining 

u = Vcp  + (V  X W)  (2) 


in  terms  of  scalar  and  vector  potentials  o>  and  W,  then  the  following  wave  equa- 
tions in  <p  and  W are  obtained  from  (1) : 


with 


= & 

Vf^W  = 

~ 3ts 


W = 0 


where 


(3) 

(4) 

(5) 


in  which  Vs  and  V0  are  the  transverse  and  dilatational  wave  velocities. 


In  our  delay  line  device,  the  transducers  are  polarized  to  provide  dis- 
placements with  a component  in  the  y-direction  only,  hence  we  assume  ux  = uz  = 0. 
Thus  the  wave  had  displacement  uy  = Uy(x,0,z,t)  in  the  y-direction  only;  the 
displacement  is  assumed  to  be  the  same  everywhere  in  y.  Thus  Eqs.(3)  and  (4) 
reduce  to: 


M. 


Vuv 


(6) 


and 


* 

4 


^2uy 

w 


az"5  > 


(7) 


Equation  (7)  is  a shear  wave  equation  with  complex  shear  modulus  p*  (one 
of  Lame's  constants).  By  the  application  of  the  method  of  separation  of 
variables  to  Eq.(7),  namely  uy(x,0,z,t)  = uy  (x,  0.  z)  e1UJt , one  obtains  for  the 
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total  solution  of  Eq.(7) 

Uy  = Ae1Ult(Ce"z+  D~BZ  ) (Eenx  +Fe-nx)  (8) 


where  uj  is  the  angular  frequency,  A,C,D,E,F  are  constants  to  be  determined  by 
the  boundary  conditions,  and  m and  n are  the  complex  propagation  constants  in 
the  z-  and  x-directions,  respectively. 


[ 


l 

L 1 

8* 


i 


Displacement  Boundary  Conditions 

This  same  equation  can  now  be  applied  to  each  of  the  three  regions  of  in- 
terest: the  solid  delay  line  and  the  upper  and  lower  liquid  layers;  then  the 
boundary  conditions  are  matched.  In  steady  state  the  frequencies  are  the  same 
in  all  regions. 

Consider  the  transverse  wave  in  the  viscoelastic  medium  above  and  below  the 
plates.  This  wave  is  traveling  unbounded  in  the  positive  z-  and  x-directions 
and  since  it  is  traveling  in  an  attenuating  medium,  the  displacement  must  approach 
zero  as  z and  x approach  infinity.  Then  C=E  = 0 in  the  liquid  above  the  plate 
and  similarly  C=F=0  in  the  liquid  below  the  plate. 

In  the  plate  itself  the  waves  are  bounded  in  the  x-direction  and  thus  there 
is  only  one  boundary  condition  that  is  C=0.  The  solutions  can  be  obtained  as 
fol lows : 


tt  —ID  z —n  x 1-CUt 

u = ° e no*  e 

Jo 

(9a) 

uyi  = H1e-“ize1U)t(enix  +e'ni*  ) 

(9b) 

u = H0e-”ozen°xe1UJt 
y 2 u 

(9c) 

where  the  subscripts  0,  1,  2 designate  the  region  above  the  plate,  the  region 
within  the  plate  and  the  region  below  the  plate,  respectively.  The  symbols  1^, 
and  Hj  are  appropriate  combinations  of  A,C,D,E,F. 


Stress  Boundary  Conditions 

In  the  presence  of  liquid  the  stress  continuity  across  the  plane  boundary, 
x=±b,  requires  that: 


U*  duv 

Ax 


(10a) 

(10b) 


1 
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If  the  plate  is  in  vacuum,  all  component  stresses  across  the  plate  surface 

Ua  AAA  Uaaaa  P AAA  „ — , , ^ H V = 4-  K AAA  AO 


must  be  zero,  hence  from  cr* 


ii  y = 0 at  x=±b  one  obtains 


e"lb  = e_nib  . 


Thus  real  [nx]  = 0 and 


IMfn,  j = ^ , 1=0,  1,  2,  .... 


If  we  now  restrict  our  consideration  to  the  solution  with  1=0  (the 
Meitzler  mode),  Eq.  (9)  for  the  wave  in  the  plate  becomes 

uv  = 2H,e-"ivzelu)t  (13) 

where  rrij  y and  m^  ^ will  refer  to  the  propagation  vector  in  the  solid,  in  the  z- 
direction,  in  the  absence  and  presence  of  liquid. 

If  the  plate  is  in  a viscoelastic  liquid,  the  complex  shear  modulus  of  the 
liquid  Ml*  is  obtained  from  the  stress  continuity  equation  (10b)  combined  with 
Eq . (9b)  to  yield : 


Mo  _ M-i 


since  (<Txy)0  = (axy)i,  at  x = b. 


% (1  - e3nib  ) 
no(l  + e2"ih) 


Next,  the  following  wave  vector  relation  for  the  liquid  is  obtained  from 
Eqs . (9a)  and  (7) : 


^ ^ = 0 . 

Mo 


From  displacement  continuity  across  the  plate,  uy  = uy^  at  x=  + b for  all 
values  of  z.  Thus  by  equating  Eqs.  (9a)  and  (9b),  one  concludes  that 


Hence  Eq.(15)  can  be  combined  with  (16)  to  obtain  an  equation  that  couples  wave 
vectors  in  the  liquid  and  in  the  solid: 

-L  -2  i UJ8  On  _ fl  l'17\ 


"gl  + "o 


"6 


Next,  the  following  wave  vector  relation  for  the  plate  (in  liquid  contact) 
is  obtained  from  Eqs. (9b)  and  (7) 


r 


n£L  + nf  + !££X  = 0 . (18) 

Hi 

While  for  the  plate  in  vacuum,  we  have  shown  by  Eq.(12)  that  n:  = 0 for  the 
Meitzler  mode,  thus  the  wave  vector  relation  (’8)  becomes 

n?¥  + ^ = 0 . (19) 

Hi 

In  our  experiment  mt  is  measured  in  air  (i.e.,  vacuum)  and  again  in  liquid. 
Thus  by  combining  Eqs.(18)  and  (19)  we  obtain 

rofi  ' “t,  = * nf  . (20) 

The  complex  modulus  of  the  viscoelastic  liquid  can  be  obtained  by  using 
the  experimental  values  of  Eq.(20)  together  with  F,qs.(l4)  and  (17)  to  give  a 
quadratic  equation  for  p,0’ . Alternatively,  an  explicit  relation  to  the  mechanical 
impedance  of  the  liquid  Z *,  can  be  obtained  as  seen  below. 

Notice  some  assumptions  that  were  made  in  the  preceding  discussion: 

(1)  Motions  are  steady-state  sinusoidal  of  small  amplitude  so  that  nonlinear 
response  of  the  propagating  media  can  be  neglected;  (2)  the  plate  is  assumed  to 
be  isotropic  and  homogeneous  and  (3)  the  transducer  will  only  excite  the  non- 
dispersive  SS(0)  mode. 

Relation  of  Mechanical  Impedance  of  Liquid  to  Propagation  Vector  in  Solid 

The  equations  derived  in  the  preceding  section  can  also  be  related  to  the 
mechanical  impedance  Z0  of  the  liquid. 

The  mechanical  impedance  of  a material  is  defined  as  the  ratio  of  the  stress 
to  the  particle  velocity  produced  by  that  stress  at  any  surface  through  the 
material.  For  the  shear  stress  components  of  interest  here  in  the  liquid  we 
can  write 


“n.  = Z Uy  n 


It  can  be  shown  that  the  characteristic  (shear)  mechanical  impedance  tor 
an  unbounded  liquid  is  related  to  the  complex  shear  modulus  and  density  by 
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Zo  = (o0UoV  . (22) 

Across  the  liquid-plate  interface,  x = + b,  we  have  seen  that  displacement 
and  stress  must  be  continuous,  hence  velocity  and  stress/velocity  must  also  be 
continuous . 

At  x = b 


°x  y 1 °xyo  Zouyo 


ax  yl  — giy°  — 2* 
“yl 


Since  we  had  for  displacement  and  stress  in  the  plate 

uyl  = 21^  cosh  n1xe_”  iz  eia)t  (24) 

CTxyi  = 21^^  nx  sinh  njXe'*  iz  eiu)t  , (25) 

the  values  of  these  functions  at  x = b are  substituted  into  Eq.(23d)  to  obtain 

zo  = - 1 tanh  ^b  . (26) 

If  we  expand  tanh  t^b,  and  take  | nx b j « 1, 

tanh  njb  = ntb  + ^Dl  L. — + ....  (27) 


z*  = i 


We  had  Eq.(20),  that  is, 


"ni  *=  L " “iv  • 


It  is  convenient  to  write  the  complex  propagation  constants  in  their  com- 
ponent form: 

mi  v = Ao  (29a) 
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mx  L = A + iB  (29b) 

where  Aq  , B0  and  A,  B are  the  attenuation  and  propagation  constants  in  the  plate 
in  the  absence  and  presence  of  liquid,  respectively.  Hence,  from  Eq. (20): 

-nf  = (A  + iB)  - (Aq  + iB0  )2  . (30) 


For  the  usual  case  where 


Eq.(30)  becomes 


where 


(B+B0)  » (A  + A0) 


-nf  = 2B0  (-AB  + iAA) 


AA  = A - Aq  and  AB  = B - B0 


and  from  Eqs.(19)  and  (29a) 


■»  * ■ 


Now  Eq.(28)  is  rewritten  as: 

Z*  = 2VS  Pj  b (AA  + iAB)  (35) 

This  is  the  basic  experimental  equation  that  relates  the  mechanical  imped- 
ance of  the  liquid  to  the  measured  differences  in  attenuation  and  phase  in  the 
plate  delay  line. 


Assume  B0  B,  and  from  Eqs.(31), 

(A  + iB)s  - (Ap  + iB0  )3  = A2  - A2  - B2  + B2  + 2i  (AB  - A0B0  ) 

= (A  -A0)(A+A0)  - (B  -B0)(B+B0)  +2i(AB  -\B0) 


= (B  + B0)(AA  A+AP.  - AB  + 2i  AA  -Jo) 

U\  R 4-  R R + R ' 


= 2B0(-AB  + iAA) 


y 

* 


I* 


r 


i 


s' 


It  is  convenient  to  summarize  the  interrelation  between  the  mechanical  im- 
pedance and  the  complex  shear  modulus  parameters  for  description  of  the  liquid. 
These  relationships  are  based  on  the  linear  theory  of  viscoelasticity.  Also 
note  that  in  viscoelastic  convention,  the  symbol  G * = p,*  for  the  shear  modulus: 


Z*  s R + iX  (36) 

G*  = g' + iG"  (37) 

G|  s (G,s  +G"2)^  (38) 

z*  = (p0G*)^  (39) 

R = pj(|G|  -G')*  (40) 

X = p|(|G|  -G')2  (41) 

G'  = (R2  -X2)/p0  (42) 

G"  = 2RX/p0  (43) 

G*  = iu)T|*  (44) 

T]*  = T]'  - iTl"  (45) 


where  R and  X are  mechanical  resistance  and  reactance,  o0  the  liquid  density, 
G'  and  G " are  the  shear  storage  and  loss  moduli,  T]'  and  T|*  are  the  real  and 
imaginary  components  of  dynamic  viscosity. 


D.  CAPILLARY  TUBE  VISCOMETRY  AND  PICNOMETRY 

The  steady-state  viscosity  and  density  of  the  test  fluids  at  atmospheric 
pressure  is  required  for  analysis  of  the  ultrasonic  data.  The  following 
tables  summarize  the  measured  values.  The  corresponding  ultrasonic  data  are 
to  be  measured  shortly. 


7 


\ 


i 


22 


fluid*-  ^'"*‘^2918) 

Oi-2-etWlfdW1  aZ'Ute 

CW°Sltl°”i  ««oSp»=« 

pressure . 


Fluid:  ATL-5067 

(Stauffer  704) 

Composition:  Tri-methylol  propane  triheptanoate 

Pressure:  1 Atmosphere 


Temperature  °C 

Viscosity  (Poises) 

Density  g/cc 

25° 

.239 

.958 

o 

CM 

.286 

.959 

15° 

.355 

.962 

10° 

.449 

.963 

SECTION  III 


CONCLUSIONS 


We  have  been  proceeding  rather  cautiously  in  trying  to  achieve  both  the 
short  and  long  term  objectives,  and  trying  to  test  each  step  along  the  way. 

But  we  do  feel  that  the  development  of  the  Diamond  Cell  Light  Scattering 
Apparatus  will  provide  an  optimum  and  unique  facility  whereby  extensive  data 
about  the  phenomenological  and  molecular  properties  of  the  lubricant  can  be 
obtained  under  laboratory  conditions  that  come  closest  to  stimulation  of  actual 
EHD  bearing  lubrication. 
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APPENDIX  A 

A NOVEL  METHOD  FOR  VISCOSITY  MEASUREMENTS  AT  HIGH  PRESSURES 


Meyer  H.  Birnboim  and  H.  Weiss 
Mechanics  Division 
Rensselaer  Polytechnic  Institute 
Troy,  New  York  12181 
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A NOVEL  METHOD  FOR  VISCOSITY  MEASUREMENTS  AT  HIGH  PRESSURES 

The  technique  of  quasi-elastic  light  scattering  is  utilized  as  a method 
to  determine  the  viscosity  of  fluids  subjected  to  large  hydrostatic  pressures. 
In  this  method,  a small  concentration  of  insoluble  polystyrene  spheres  of  known 
diameter,  0.109  microns,  were  added  to  the  fluid,  water.  The  high  pressure 
cell  in  which  the  fluid  is  contained  has  three  optical  windows.  The  laser 
light  scattered  from  the  suspended  spheres  is  analyzed  to  determine  the  diffu- 
sion coefficient  of  the  spheres.  Since,  by  Stokes  law, 

the  viscosity  of  the  fluid  is  determined  from  the  measured  diffusion  coeffi- 
cient and  known  particle  radius  r,  A small  correction  is  required  for  the 
contraction  of  the  sphere  and  for  refractive  index  change  of  the  water  under 
pressure.  Results  on  the  viscosity  of  water  at  25°C  up  to  pressures  of 
2.5  kilobars  is  reported  and  compared  to  published  results. 

The  high  pressure  system  and  scattering  cell  (1)  is  illustrated  in 
Figure  1.  The  system  is  assembled  from  standard  high  pressure  components. 

The  scattering  cell  consists  of  a high  pressure  cross;  three  blank  conical 
plugs  are  drilled  out  and  fitted  with  sapphire  cylinder  windows  which  are 
epoxied  in  place.  Windows  from  1/8"  to  1/4"  diameter  by  3/4"  long  have  been 
used  successfully  to  pressures  of  35,000  psi.  The  fluid  is  filled  through  the 
reservoir  to  overflow  at  the  window  to  ensure  that  no  air  is  trapped,  then 
the  window  plug  is  tightened.  Pressure  is  generated  by  hand,  by  compression 


American  Instrument  Company,  Silver  Springs,  Maryland. 
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of  fluid  in  the  valves.  The  sequence  is:  close  all  valves,  open  valves  A,  B, 

close  A,  open  C,  close  B,  C,  then  repeat  sequence  from  open  valves  A,  B. 

Three  repetitions  of  the  sequence  will  bring  the  pressure  to  35,000  psi  for 
a water  sample. 

The  light  scattered  by  the  polystyrene  spheres  is  observed  at  a scattering 
angle  of  90°  to  the  main  beam.  The  intensity  of  scattered  light  is  monitored 
by  a photomultiplier  which  is  operated  in  a photon  counting  mode,  thus  the 
number  of  photoelectric  events  (photocounts)  per  second  n(t)  is  proportional 
to  the  intensity  of  light.  The  intensity  is  not  constant  on  short  time  scales 
but  exhibits  fluctuation  associated  with  the  Brownian  motion  of  the  scattering 
particles  [see  review  by  Cummins  and  Swinney  (2)].  The  sequence  of  "instantane- 
ous" intensities  are  stored  in  real  time  in  the  memory  of  a computer,  then  the 
autocorrelation  function  is  calculated. 

The  photocount  rate  autocorrelation  function  c(t)  is  defined  in  terms  of 
the  sequence  n(t) ; and  in  turn  c(t)  is  related  (2)  to  the  translational  diffu- 
sion coefficient  D: 

2 

0 - ?TW 

c(t)  = <n(t)n(t  +T))t  = <n(t)>  6(t)  + <n(t)>t  + Ae  T (2) 


The  brackets  ( ) denote  the  average  overall  starting  times  t;  6(t)  = 1 for 
t = 0 and  5(t)  = 0 otherwise. 

For  our  computer  method,  the  real  time  t and  delay  time  t are  always  a 
multiple  of  T,  that  is  t = jT  with  j = 0,1,2  ....  ; A is  a coefficient  related 
to  the  intensity  of  light  scattered  from  the  macromolecules;  and  the  magnitude 
of  the  scattering  vector  K is  obtained  from 


K 


2 


4Tm  . 2 

■17*) 

o 


(3) 
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where  n and  \ are  the  refractive  index  of  the  solution  and  wavelength  of  the 
o o 

light  in  vacuum.  In  this  experiment,  the  scattering  angle  9 was  kept  at  90° . 

2 

The  correlation  time,  1/2  DK  , and  therefore  D,  is  calculated  from  the  slope 
of  inf" c (t)  - (n(t))2l  vs.  r,  where  the  shot  noise,  the  first  term  of  Eq.(2) 
may  be  neglected  since  it  is  zero  except  at  t = 0,  and  (n(t))t  is  the  average 
count  rate . 

An  example  of  the  autocorrelation  function  is  seen  in  Figure  2,  taken  at 
a pressure  of  15,000  psi.  The  diffusion  coefficient  determined  from  the  slope 
of  this  curve  in  conjunction  with  Eq.(l)  yields  the  viscosity  at  15,000  psi. 
Our  results  for  the  limited  pressure  range  of  0 to  35,000  psi  (2.5  kilobars) 
taken  at  a temperature  of  23°C  are  shown  by  the  points  in  Fig. 3.  These  points 
are  superposed  on  the  extensive  data  of  Butt  and  Capi  (3)  as  represented  by 
the  solid  lines.  The  results  are  only  intended  as  an  illustrative  example  of 
the  method. 
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FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Block  diagram  of  the  Homodyne  Spectrometer  and  high  pressure 
cell  for  the  quasi-elastic  light  scattering  experiment. 

Light  scattered  at  an  angle  of  90°  is  detected  by  the  photo- 
multiplier. The  number  of  photoelastic  pulses  n^  in  each 
successive  interval  of  length  T are  stored  in  the  memory  of 
the  computer. 


Plot  of  fe  I <n.n.,.>.  - <n.>.  I vs.  delay  time  t . at  a 
l i+j  i li  7 j 


scattering  angle  of  90°  and  atmospheric  pressure.  The  photo- 


pulse autocorrelation  function  <n^n^+^ > ^ and  the  average 


number  of  counts  are  calculated  from  the  measured  sequence 
n..  The  t.  = jT,  where  j = 0,1,2,...  . As  discussed  in  the 

-1  _ 


J 


text,  the  slope  of  the  above  curve  yields  t = 2D  K , where 
’ r c tt  ’ 


t is  the  characteristic  correlation  time, 
c 


Viscosity  of  water  vs.  pressure.  The  filled  points  and  solid 
lines  are  data  of  Bett  and  Cappi.  The  open  circles  and 
dashed  lines  are  our  data. 


Figure  1 Block  Diagram  of  the  Homodyne  Spectrometer. 
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Figure  3 Viscosity  of  Water  vs.  Pressure 
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It  has  been  shown  experimentally  by  various  investigators  that  the 
dynamic  viscosity  of  a dilute  polymer  solution  in  the  limit  of  high  fre- 
quency does  not  become  the  viscosity  of  the  solvent.  To  study  this 
limiting  behavior,  a high  frequency  surface  delay  line  has  been  developed 
for  measurement  of  and  in  the  range  from  2 MHz  to  40  MHz.  Details 
and  equations  for  this  device  are  described  herein.  Also  the  special  data 
acquisition  technique  required  to  extract  phase/amplitude  information  from 
RF  pulses  is  described. 

Three  Newtonian  fluids,  NBS  oil  D,  NBS  oil  J and  NBS  oil  L are  used 
to  obtain  calibration  curves  which  confirm  the  theoretical  analysis.  In 
addition  data  is  presented  on  one  viscoelastic  fluid,  NBS  oil  N. 


1 . Introduction 


The  viscoelastic  response  of  dilute  polymer  solutions  is  usually  separated 
into  two  parts,  the  polymer  and  the  solvent  contributions.  The  theories  of 
Rouse,  Zimm  and  others  attribute  the  polymer  contribution  to  internal  relaxa- 
tional  modes  of  the  macromolecule,  the  contribution  to  the  total  viscosity  at 
low  frequencies  drops  to  zero  in  the  limit  of  high  frequencies.  On  this  basis, 
one  expects  the  viscosity  of  the  solution  to  approach  the  viscosity  of  the 
solvent  in  the  limit  of  high  frequencies. 

Several  investigators  have  found  that  this  assumption  is  not  borne  out 
experimentally,  and  the  difference  between  the  limiting  solution  viscosity  and 
that  of  the  solvent  has  been  attributed  to  internal  friction.  This  inter- 
pretation has  recently  been  brought  into  question  by  Schrag  and  the  problem 
remains  open. 

The  experimental  techniques  that  have  been  used  to  investigate  this  phe- 
nomenon have  primarily  been  in  the  frequency  range  from  0.01  Hz  to  10  KHz. 

This  means  that  it  has  been  necessary  to  use  highly  viscous  solvents  together 
with  the  principle  of  time-temperature  superposition  to  obtain  the  limiting 
behavior.  Wide  frequency  range  measurements  at  a single  temperature  would  be 
desirable  to  obtain  a less  ambiguous  interpretation  of  results. 

In  this  paper  we  describe  a device  and  method  for  measuring  the  visco- 
elastic response  of  low  viscosity  liquids  in  the  frequency  range  2 MHz  to 
40  MHz. 

The  device  consists  of  a glass  plate  in  which  a shear  wave  is  propagated 
between  the  ends,  while  the  liquid  is  in  contact  with  the  major  surfaces.  The 
transverse  mode  of  vibration,  called  the  Meitzler  mode,  is  the  same  mode  as  was 
used  by  Meeker  and  Meitzler  in  an  aluminum  strip  delay  line,  and  was  subsequently 
applied  by  Hunston  and  Meyers  to  polymer  solutions.  The  technique  here  differs 
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from  the  high  frequency  technique  of  McSkimin  in  that  the  wave  propagates  in 
the  plane  of  the  solid/liquid  interface  instead  of  at  a normal  or  oblique 
incidence  to  the  interface  as  in  the  two  methods  of  McSkimin. 

2 . Description  of  Apparatus 
a)  The  Surface  Delay  Line 

The  geometry  of  the  fused  silica  surface  delay  line  is  seen  in  Fig.l. 

A transverse  wave  (Meitzler  mode)  polarized  in  the  y-direction  is  propagated 
in  the  z-direction.  Piezoelectric  ceramic  driving  and  response  transducers 
of  center  frequency  25  MHz,  are  fixed  at  z = 0 and  z - l,  respectively.  The 
line  length  & = 7.62  cm  and  thickness  2b  = 0.158  cm.  The  liquid  is  in  contact 
with  the  major  surfaces  at  the  interface  x = + b.  The  beam  (transducer) 
width  is  1.9  cms , and  the  width  of  the  line  is  5.08  cm.  The  thickness  of 
liquid  layer  for  the  data  reported  here  is  0.025  cm.  The  line  is  isolated 
from  its  polycarbonate  case  by  0-rings.  and  has  a flow  through  arrangement  for 
filling . 

* 

The  coupled  wave  solution  for  displacement  in  the  y-direction  u,  in  the 
solid  and  in  the  liquid  are  respectively  given  by: 

I / 2 v£  1 


i : 

( 2 

W J 

ug(x,z,t)  = Hg  exp  -i(^2  " 

i “*  ,)  z + iout  | 

f.  1 

bDV  J 

r .(m  . uuZ*  1 f . ujpL  ] 

- -i  ^ - i i z + lout  I exp;  - i — x ! 


uL(x.z,t)  = exp  - ! - i^  2 * i ^2/  z + ioutj  exP[ 


Through  Eq.  (1).  the  difference  in  attenuation  £A  and  in  phase  AB  of  the  wave 

measured  at  z - jt  in  the  absence  and  in  the  presence  of  liquid  can  be  related 

¥ 

to  the  mechanical  impedance  Z of  the  liquid  by: 

¥ 

Complete  equations  developed  in  the  Appendix  are  not  included  here. 
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(3) 


M + iAB  = 2^b?  Z< 

where  o andc;  are  the  density  of  the  solid  and  the  liquid  and  V is  the 

velocity  in  the  solid  without  liquid. 

* _ 

The  mechanical  impedance  Z = R + iX  is  related  to  the  dynamic  shear 
moduli  G = G'  + iG"  = G'  + icuT]'  by: 

OjG'  = R2  - X2  ; plG"  = 2RX  (4) 

For  a Newtonian  liquid  of  viscosity  T|,  G'  = 0,  thus 


b)  Pulse  Technique  and  Data  Acquisition 

A cross-correlation  method  for  measuring  phase  shift  and  amplitudes 
of  radio  frequency  (RF)  continuous  wave  signals,  called  the  methof  of  walking 
frequencies,  has  been  described  previously  (1).  Herein  this  method  is  ex- 
tended to  phase/amplitude  measurements  on  pulsed  RF  echoes. 

As  is  seen  in  Fig. 2,  a 1 MHz  crystal  clock  is  used:  to  drive  the  syn- 

thesizer to  produce  the  frequency  of  f MHz,  to  drive  a programmable  clock  to 
produce  the  pulse  P every  40  ps,  and  to  drive  a snap  diode  which  on  mixing 
with  P produces  a very  fast  trigger  pulse  T every  40  ps  . The  RF  drive  pulse 
A of  frequency  f is  produced  every  1280  ps  (32  X 40  ps)  by  mixing  the  synthe- 
sizer output  with  the  d.c.  pulse  Q of  nominal  width  25  ps , that  has  been 
appropriately  delayed.  The  response  signals  (echoes)  B,  from  the  delay  line 
output  are  then  sampled  every  40  ps  by  the  trigger  pulse  T.  The  length  of 
line  is  such  that  the  interval  between  successive  echoes  is  nominally  40  ps ; 
thus  the  drive  delay  (above)  is  adjusted  so  that  echoes  and  sampling  overlap. 

The  output  of  the  sample  and  hold  amplifier  is  then  measured  by  the 
analog  to  digital  converter  every  40  ps  and  stored  into  the  memory  of  the 
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4. 


computer.  In  order  that  a complete  and  separate  sinusoidal  wave  is  developed 
in  the  memory  from  each  echo,  the  synthesizer  must  always  be  set  to  a walking 
frequency.  A walking  frequency  f is  given  by: 

N 1 

* fN  + t ~ T + 32 nT  (6) 

where  f^  corresponds  to  any  commensurate  frequency  (every  25  KHz  since 
T = 40  ns  and  N = 1,2,3  ...)  and  if  we  choose  n = 625  memory  locations  to 
correspond  to  2tt  radians,  then  Af  = 1.25  Hz.  Thereby  for  each  echo  one 
sinusoidal  wave  is  reconstructed  in  625  locations,  and,  if  they  existed, 

32  echoes  could  be  measured. 

By  the  cross-correlation  techniques  described  elsewhere  (2)  the  relative 
phase  and  amplitude  between  any  one  echo  and  another  echo  or  drive  signal 
can  be  computed. 

3.  Experimental  Results  and  Conclusions 

As  can  be  seen  from  Eq.  (3)  we  expect  linear  relations  between  AA  and  R 
and  AB  and  X,  with  the  proportionality  constant  indicated.  Alternatively 
calibration  curves  can  be  obtained  by  using  Newtonian  fluids  for  which  R 
and  X can  be  calculated  by  Eq.  (5),  from  steady-state  viscosity,  frequency 
and  density.  Such  calibration  curves  are  shown  in  Fig. 3.  The  superposition 
of  these  various  liquids  and  frequencies  establish  the  validity  of  the  method. 

One  fluid,  NBS  oil  N,  is  shown  in  Fig. 4 to  be  viscoelastic  in  this  fre- 
quency range. 
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FIGURE  CAPTIONS 


Figure  1 


Figure  2 


Figure  3 


Figure  4 


Delay  line  configuration,  where  k is  the  propagation 

vector,  u is  the  particle  displacement,  A and  B are  the 
drive  and  response  RF  pulses. 

Block  diagram  of  electronics.  The  output  of  the  analog 
to  digital  converter  CPU,  transfers  data  directly  into 
computer  memory. 

Calibration  curves  for  liquid  mechanical  resistance  and 
reactance  in  terms  of  measured  attenuation  AA  and  phase 
shift  AB,  respectively.  Measurements  at  3.4,  6.0  and 
14.5  MHz. 

Viscoelastic  behavior  of  NBS  oil  N. 
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Figure  1 Delay  Line  Configuration  where  k is  the  Propagation  Vector, 
u is  the  Particle  Displacement,  A and  B are  the  Drive  and 
Pesponae  RF  Pulses 


CPU 


Figure  2 Block  Diagram  of  Electronics.  The  output  of  the  analog  to  digital 
converter  CPU,  transfers  data  directly  Into  computer  memory. 
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Figure  1(a) 


Stress/Def ormation  History  in  the  EHD  Problem 
(Contact  Zone) 
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Figure  1(b)  Stress/Deformation  History  in  the  EHD  Problem 
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Figure  2 Survey  of  Properties  of  the  Liquid  State 
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Figure  3 The  Volume-Temperature  Curves  with  Reference  to: 

Equilibrium.  Nonequilibrium.  Glass/I. iquid  to  Glass 
T rans i t ion 
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Figure  4(b)  The  Dynamic  Shear  Viscosity  at  Thermodynamic 
Equilibrium  and  Atmospheric  Pressure 
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Figure  5(a)  The  Dynamic  Compressibility:  Viscosity 

and  Modulus 
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Figure  5(b)  The  Dynamic  Compressibility:  Viscosity 
and  Modulus 
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Figure  9 The  New  High  Pressure  Cell 
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Figure  12  - Captions 

A - Argon  ion  laser,  B - Single  frequency  etalon,  C - Polarization  rotator, 

D - Spatial  filter,  F.  - Partial  (5%)  reflector  - for  beam  intensity  monitor 
or  alternately  Es  entrance  point  for  HeNe  (red)  laser  (second  source), 

F - Lens  to  focus  beam  into  "waist"  diameter  of  50  microns,  G - Prism.  Rotate 
180°  to  displace  beam  from  "forward  to  backward  scattering"  illumination, 

H - Diamond  anvil  cell,  I - Cell  clamp  with  cell  orientation  adjustment, 

J - Pressure  adjustment  screw  (only  partly  shown),  K - Pneumatic  impact  relay, 
(not  shown),  L - Primary7  goniometer  to  adjust  scattering  angle,  M,  - Mirror 
for  beam  steering  in  forward  scattering  regime,  H - Partially  reflecting  mirror 
for  beam  steering  in  back  scattering  regime,  N - Adjustable  slits  or  apertures, 

0 - Microscope  objective,  P - Three  port  microscope  body  attached  to  Table  L, 

M,^  - Rotatable  and  removable  mirror  for  beam  steering  to  R,  S;  R - Eyepiece 
or  Vidicon  port,  S - Lens,  T - Piezoelectrically  driven  Fabry  Perot  interfero- 
meter of  adjustable  free  spectral  range,  U - Multipass  corner  cube  reflector 
and  an  Iodine  filter  (not  indicated),  V - Photomultiplier.  W - Lens,  X - Double 
grating  spectrometer,  Y - Filter  for  fluorescence. 
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